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Abstract 

In the Standard Model of particle physics, the strength of the couplings of the b quark to 
the u and c quarks, \V u b\ and \V c b\, are governed by the coupling of the quarks to the Higgs 
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for the A® baryon to decay into the pp~V^ final state relative to the A+final state is 
measured. Combined with theoretical calculations of the strong interaction and a previously 
measured value of |U c fe|) the first \V u b\ measurement to use a baryonic decay is performed. 
This measurement is consistent with previous determinations of \V u t,\ using B meson decays 
to specific final states and confirms the existing incompatibility with those using an inclusive 
sample of final states. 
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uncertainty is dominated by the LQCD calcula¬ 
tions, which have recently been updated 12 13 


In the Standard Model (SM) of particle physics, 
the decay of one quark to another by the emis¬ 
sion of a virtual W boson is described by the 3x3 
unitary Cabibbo-Kobayashi-Maskawa (CKM) 
matrix 11,21. This matrix arises from the cou¬ 
pling of the quarks to the Higgs boson. While 
the SM does not predict the values of the four 
free parameters of the CKM matrix, the mea¬ 
surements of these parameters in different pro¬ 
cesses should be consistent with each other. If 
they are not, it is a sign of physics beyond the 
SM. In global fits combining all available mea¬ 
surements |:3||4j, the sensitivity of the overall 
consistency check is limited by the precision in 
the measurements of the magnitude and phase 
of the matrix element V u b, which describes the 
transition of a 6 quark to a u quark. 

The magnitude of 14& can be measured 
via the semileptonic quark-level transition 
6 —> ut~Vt. Semileptonic decays are used to 
minimise the uncertainties arising from the in¬ 
teraction of the strong force, described by quan¬ 
tum chromodynamics (QCD), between the final- 
state quarks. For the measurement of the mag¬ 
nitude of Vub, as opposed to measurements of 
the phase, all decays of the b quark, and the 
equivalent 6 quark, can be considered together. 

There are two complementary methods to per¬ 
form the measurement. From an experimental 
point of view, the simplest is to measure the 
branching fraction (probability to decay to a 
given final state) of a specific (exclusive) decay. 

An example is the decay of a B° ( bd ) meson 
to the final state tt + £~u, where the influence 
of the strong interaction on the decay, encom¬ 
passed by a B° —> 7r + form factor, is predicted | Vhb | 2 /1 | 2 according to 


and result in larger values of V u b than the average 
given in Ref. |7|. The alternative method is to 
measure the differential decay rate in an inclusive 
way over all possible B meson decays contain¬ 
ing the b —>- ui~i 7 quark level transition. This 
results in |I4&| = (4.41 ±0.15l°u7) x 10 -3 |14], 
where the first uncertainty arises from the ex¬ 
perimental measurement and the second from 
theoretical calculations. The discrepancy be¬ 
tween the exclusive and inclusive 14,6 determi¬ 
nations is approximately three standard devi¬ 
ations and has been a long-standing puzzle in 
flavour physics. Several explanations have been 
proposed, such as the presence of a right-handed 
(vector plus axial-vector) coupling as an exten¬ 
sion of the SM beyond the left-handed (vector 
minus axial-vector) W coupling 15-18J. A simi¬ 


lar discrepancy also exists between exclusive and 
inclusive measurements of 14 , (the coupling of 
the b quark to the c quark) [fT. 

This article describes a measurement of the 
ratio of branching fractions of the (bud) 
baryon into the p£~ v and A+£~v final states. 
This is performed using proton-proton collision 
data from the LHCb detector, corresponding to 
2.0 fb -1 of integrated luminosity collected at a 
centre-of-mass energy of 8 TeV. The b —> u tran¬ 
sition, T° —>• pp~u^, has not been considered be¬ 
fore as A ° baryons are not produced at an e + e“ 
R-factory; however, at the LHC, they consti¬ 
tute around 20% of the 6-hadrons produced [l9|. 
These measurements together with recent LQCD 
calculations 1201 allow for the determination of 


by non-perturbative techniques such as lattice 
QCD (LQCD) |5| or QCD sum rules |6|. The 
world average from Ref. |7| for this method, us¬ 
ing the decays B° —>■ 7r + £~u and B~ —> tt°£~u, 
is \V u b\ = (3.28 ± 0.29) x 10~ 3 , where the 
most precise experimental inputs come from the 
BaBar |8.9 and Belle [10,11 experiments. The 
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where B denotes the branching fraction and Rp p 
is a ratio of the relevant form factors, calcu¬ 
lated using LQCD. This is then converted into a 
measurement of \V u b\ using the existing measure- 
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merits of |V^| obtained from exclusive decays. 
The normalisation to the A® —> T+ li~V^ decay 
cancels many experimental uncertainties, includ¬ 
ing the uncertainty on the total production rate 
of T° baryons. At the LHC, the number of signal 
candidates is large, allowing the optimisation of 
the event selection and the analysis approach to 
minimise systematic effects. 


The LHCb detector 21 22 is one of the four 
major detectors at the Large Hadron Collider. 
It is instrumented in a cone around the proton 
beam axis, covering the angles between 10 and 
250 mrad, where most 6-hadron decays produced 
in proton-proton collisions occur. The detector 
includes a high-precision tracking system with 
a dipole magnet, providing a measurement of 
momentum and impact parameter (IP), defined 
for charged particles as the minimum distance of 
a track to a primary proton-proton interaction 
vertex (PV). Different types of charged particles 
are distinguished using information from two 
ring-imaging Cherenkov detectors, a calorime¬ 
ter and a muon system. Simulated samples of 
specific signal and background decay modes of 
6 hadrons are used at many stages throughout 
the analysis. These simulated events model the 
experimental conditions in full detail, including 
the proton-proton collision, the decay of the par¬ 
ticles, and the response of the detector. The 
software used is described in Refs. 


23-29 


Candidates of the signal modes are required 
to pass a trigger system 1301 which reduces in 
real-time the rate of recorded collisions (events) 
from the 40 MHz read-out clock of the LHC to 
around 4 kHz. For this analysis, the trigger re¬ 
quires a muon with a large momentum transverse 
to the beam axis that at the same time forms 
a good vertex with another track in the event. 
This vertex should be displaced from the PVs in 
the event. The identification efficiency for these 
high momentum muons is 98%. 

In the selection of the final states, stringent 
particle identification (PID) requirements are 


applied to the proton. These criteria are accom¬ 
panied by a requirement that its momentum is 
greater than 15 GeV/c as the PID performance is 
most effective for protons above the momentum 
threshold to produce Cherenkov light. The ppL~ 
vertex fit is required to be of good quality, which 
reduces background from most of the 6 —s- cp,~v^ 
decays as the resulting ground state charmed 
hadrons have significant lifetime. 

To reconstruct —> (T+ — >■ pK~7T + )fi~u fl 
candidates, two additional tracks, positively 
identified as a pion and kaon, are combined 
with the proton to form a A+ —)■ pK~n + can¬ 
didate. These are reconstructed from the same 
pp~ vertex as the A® —>■ ppru^ signal to min¬ 
imise systematic uncertainties. As the lifetime 
of the A+ is short compared to other weakly 
decaying charm hadrons, the requirement has 
an acceptable efficiency. 

There is a large background from 6-hadron 
decays with additional charged tracks in the de¬ 
cay products, as illustrated in Fig. [lj To reduce 
this background, a multivariate machine learning 
algorithm (a boosted decision tree, BDT |31,321) 
is employed to determine the compatibility of 
each track from a charged particle in the event 
to originate from the same vertex as the signal 
candidate. This isolation BDT includes vari¬ 
ables such as the change in vertex quality if 
the track is combined with the signal vertex, 
as well as kinematic and IP information of the 
track that is tested. For the BDT, the train¬ 
ing sample of well isolated tracks consists of 
all tracks apart from the signal decay products 
in a sample of simulated A° —> events. 

The training sample of non-isolated tracks con¬ 
sists of the tracks from charged particles in 
the decay products X in a sample of simulated 
T° —> (%+ —> pX)ijrv IM events. The BDT se¬ 
lection removes 90% of background with addi¬ 
tional charged particles from the signal vertex 
while it retains more than 80% of signal. The 
same isolation requirement is placed on both the 
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Figure 1: Diagram illustrating the topol¬ 
ogy for the (top) signal and (bottom) back¬ 
ground decays. The A® baryon travels about 1 cm 
on average before decaying; its flight direction is 
indicated in the diagram. In the 7l °—y pp-jr^ signal 
case, the only other particles present are typically 
reconstructed far away from the signal, which are 
shown as grey arrows. For the background from A+ 
decays, there are particles which are reconstructed 
in close proximity to the signal and which are indi¬ 
cated as dotted arrows. 


A ° —» pp~V^ and A ° —>- (%+ —>- pK~Ti + )p~u ll 
decay candidates, where the pion and kaon are 
ignored in the calculation of the BDT response 
for the yl°— > (A+ — > pK~7T + )p~u M case. 

The A® mass is reconstructed using the so- 
called corrected mass [33], dehned as 

m CO rr = \firf-+p i + p±, 

where is the visible mass of the hp pair and 
p± is the momentum of the hp pair transverse 


to the A® flight direction, where h represents 
either the proton or T+ candidate. The flight 
direction is measured using the PV and A® vertex 
positions. The uncertainties on the PV and 
the A® vertex are estimated for each candidate 
and propagated to the uncertainty on m corr ; the 
dominant contribution is from the uncertainty 
in the A® vertex. 

Candidates with an uncertainty of less than 
100MeV/c 2 on the corrected mass are selected 
for the T° —> pp~V^ decay. This selects only 23% 
of the signal; however, the separation between 
signal and background for these candidates is 
significantly improved and the selection thus re¬ 
duces the dependence on background modelling. 

The LQCD form-factors that are required to 
calculate \V ub \ are most precise in the kinematic 
region where q 2 , the invariant mass squared of 
the muon and the neutrino in the decay, is high. 
The neutrino is not reconstructed, but q 2 can 
still be determined using the A® flight direction 
and the T° mass, but only up to a two-fold 
ambiguity. The correct solution has a resolu¬ 
tion of about 1 GeV 2 /c 4 , while the wrong solu¬ 
tion has a resolution of 4GeV 2 /c 4 . To avoid 
influence on the measurement by the large un¬ 
certainty in form factors at low g 2 , both so¬ 
lutions are required to exceed 15 GeV 2 /c 4 for 
the T° —> pp~Uf, decay and 7GeV 2 /c 4 for the 
A° b ^ (A+ -+pK~ '7 r + )p decay. Simulation 

shows that only 2% of T° —> pp~v^ decays and 
5% of —> A+ /i~z% decays with q 2 values below 
the cut values pass the selection requirements. 
The effect of this can be seen in Fig. [2] where the 
efficiency for signal below 15GeV 2 /c 4 is reduced 
significantly if requirements are applied on both 
solutions. It is also possible that both solutions 
are imaginary due to the limited detector resolu¬ 
tion. Candidates of this type are rejected. The 
overall q 2 selection has an efficiency of 38% for 
T° —> pp~v^ and 39% for T° —>• kt+decays 
in their respective high-g 2 regions. 

The mass distributions of the signal candi- 
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Figure 2: Illustrating the method used to re¬ 
duce the number of selected events from the 
q 2 region where lattice QCD has high uncer¬ 
tainties. The efficiency of simulated AQ —?• pp~u ll 
candidates as a function of q 2 . For the case where 
one q 2 solution is required to be above 15 GeV 2 /c 4 
(marked by the vertical line), there is still significant 
efficiency for signal below this value, whereas, when 
both solutions have this requirement, only a small 
amount of signal below 15GeV 2 /c 4 is selected. 
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dates for the two decays are shown in Fig. [3] The 
signal yields are determined from separate y 2 
fits to the m corr distributions of the T° —> ppTv^ 
and/lj)—(T+ —>• pK~7T + )p~u fl candidates. The 
shapes of the signal and background components 
are modelled using simulation, where the un¬ 
certainties coming from the finite size of the 
simulated samples are propagated in the fits. 
The yields of all background components are 
allowed to vary within uncertainties obtained as 
described below. 

For the fit to the m corr distribution of the 
A ° —> pp~V^ candidates, many sources of back¬ 
ground are accounted for. The largest of 
these is the cross-feed from A® —s- T+ p~v tl 
decays, where the T+ decays into a pro¬ 
ton and other particles that are not recon¬ 
structed. The amount of background arising 


Figure 3: Corrected mass fit used for de¬ 
termining signal yields. Fits are made to 
(top) A® -» pp~V^, and (bottom) A ° —> (T+ —> 
pK~Ti + )p,~V^ candidates. The statistical uncer¬ 
tainties arising from the finite size of the simulation 
samples used to model the mass shapes are indi¬ 
cated by open boxes while the data are represented 
by the black points. The statistical uncertainty on 
the data points is smaller than the marker size used. 
The different signal and background components 
appear in the same order in the fits and the legends. 
There are no data above the nominal A® mass due 
to the removal of unphysical q 2 solutions. 

from these decay modes is estimated by fully 
reconstructing two T+ decays in the data. The 
background where the additional particles in¬ 
clude charged particles originating directly from 
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the A+ decay is estimated by reconstructing 
A ®—> (A+ — » pK~n + )p~T' lx decays, whereas the 
background where only neutral particles come 
directly from the A+ decay is estimated by 
reconstructing yl° —> (yl+ —* pK®)p~i decays. 
These two background categories are separated 
because the isolation BDT significantly reduces 
the charged component but has no effect on 
the neutral case. For the rest of the A+ de¬ 
cay modes, the relative branching fraction be¬ 
tween the decay and either the A+ —> pK~n + 
or A+ —> pKg decay modes, as appropriate, is 
taken from Ref. |14j. For some neutral decay 
modes, where only the corresponding mode with 
charged decay particles is measured, assump¬ 
tions based on isospin symmetry are used. In 
these decays, an uncertainty corresponding to 
100% of the branching fraction is allowed for 
in the fit. Background from A® —> D 0 pp~u fJ/ 
decays is constrained in a similar way to the 
A+ charged decay modes, with the normalisa¬ 
tion done relative to A ° — > D°(— > K~Ti + )pp~V ^ 
decays reconstructed in the data. 

Any background with a A+ baryon may 
also arise from decays of the type A b —* 
(A* + —* A^7nr)p~T' ll , where A* c + represents the 
A c (2595) + or A c (2625) + resonances as well as 
non-resonant contributions. The proportions 
between the A® — > ( A* + —>■ A^ti’k)p~V 11 and the 
^b ~‘" A+p, backgrounds are determined from 
the fit to the A b — >■ (A+ — > pK~n + )p~T' fi m co „ 
distribution and then used in the A b —> pp~u fl 
fit. 

The decays A b —>■ N*p~v^, where the N* 
baryon decays into a proton and other non- 
reconstructed particles, are very similar to the 
signal decay and have poorly known branch¬ 
ing fractions. The N* resonance represents any 
of the states 1V(1440), 1V(1520), 1V(1535) or 
IV(1770). None of the A b —> N*p~u fl decays 
have been observed and the m corr shape of these 
decays is obtained using simulation samples gen¬ 
erated according to the quark-model prediction 


of the form factors and branching fractions 34 


A 100% uncertainty is allowed for in the branch¬ 
ing fractions of these decays. 

Background where a pion or kaon is mis- 
identified as a proton originates from various 
sources and is measured by using a special data 
set where no PID is applied to the proton can¬ 
didate. Finally, an estimate of combinatorial 
background, where the proton and muon origi¬ 
nate from different decays, is obtained from a 
data set where the proton and muon have the 
same charge. The amount and shape of this 
background are in good agreement between the 
same-sign and opposite-sign pji samples for cor¬ 
rected masses above 6 GeV/c 2 . 

For the A b —>- (A+ —> pK~7r + )p~u fi yield, 
the reconstructed pK~7T + mass is studied to de¬ 
termine the level of combinatorial background. 
The A~ 1" signal shape is modelled using a Gaus¬ 
sian function with an asymmetric power-law tail, 
and the background is modelled as an exponen¬ 
tial function. Within a selected signal region of 
30MeV/c 2 from the known A+ mass the combi¬ 
natorial background is 2% of the signal yield. 
Subsequently, a fit is performed to the m corr 
distribution for A b — >■ (A+ — pK~7T + )p~u /1 
candidates, as shown in Fig. [3| which is used 


to discriminate between A b —s- A+p 


and 


A b — > ( A* + —>■ A^'K’k)p~V ^ 1 decays. 

The A b —> pp-v^ and A b —> (A+ — > 
pK~ 7r + )/i _ I% yields are 17,687±733 and 
34,255±571, respectively. This is the first obser¬ 
vation of the decay A b —>• pp~V^. 

The A b —>■ pp^u^ branching fraction is mea¬ 
sured relative to the A b — >• (A+ — >• pK~ir + )p~v M 
branching fraction. The relative efficiencies for 
reconstruction, trigger and final event selection 
are obtained from simulated events, with several 
corrections applied to improve the agreement 
between the data and the simulation. These 
correct for differences between data and simu¬ 
lation in the detector response and differences 
in the A b kinematic properties for the selected 
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A b —> pn~Vp and A b —y (A+ —>• pK~n + )p~T'^ 
candidates. The ratio of efficiencies is 3.52±0.20, 
with the sources of the uncertainty described be¬ 
low. 

Systematic uncertainties associated with the 
measurement are summarised in Table □3 The 
largest uncertainty originates from the A+ —> 
pK~n + branching fraction, which is taken from 
Ref. 35 . This is followed by the uncertainty 


on the trigger response, which is due to the 
statistical uncertainty of the calibration sam¬ 
ple. Other contributions come from the track¬ 
ing efficiency, which is due to possible differ¬ 
ences between the data and simulation in the 
probability of interactions with the material 
of the detector for the kaon and pion in the 
21°—)- (A+ —» pK~ 7r + ) n~Vn decay. Another sys¬ 
tematic uncertainty is assigned due to the lim¬ 
ited knowledge of the momentum distribution 
for the 7l+ —> pK~n + decay products. Uncer¬ 
tainties related to the background composition 
are included in the statistical uncertainty for 
the signal yield through the use of nuisance pa¬ 
rameters in the fit. The exception to this is the 
uncertainty on the A® —>■ N*p~V^ mass shapes 
due to the limited knowledge of the form factors 
and widths of each state, which is estimated by 
generating pseudoexperiments and assessing the 
impact on the signal yield. 

Smaller uncertainties are assigned for the 
following effects: the uncertainty in the A® life¬ 
time; differences in data and simulation in the 
isolation BDT response; differences in the rel¬ 
ative efficiency and q 2 migration due to form 
factor uncertainties for both signal and normali¬ 
sation channels; corrections to the A® kinematic 
properties; the disagreement in the q 2 migra¬ 
tion between data and simulation; and the finite 
size of the PID calibration samples. The to¬ 
tal fractional systematic uncertainty is Ig'®%, 
where the individual uncertainties are added in 
quadrature. The small impact of the form factor 
uncertainties means that the measured ratio of 


Tabic 1: Summary of systematic uncertainties. 

The table shows the relative systematic uncertainty 
on the ratio of the A b —> pp^V^ and A® —> 2l+p~z7 M 
branching fractions broken into its individual con¬ 
tributions. The total is obtained by adding them in 
quadrature. Uncertainties on the background levels 
are not listed here as they are incorporated into the 
fits. 


Source Relative uncertainty (%) 

—>■ pK + ir ) 

+ 4.7 

- 5.3 

Trigger 

3.2 

Tracking 

3.0 

A+ selection efficiency 

3.0 

A b -y N* p~l shapes 

2.3 

A b lifetime 

1.5 

Isolation 

1.4 

Form factor 

1.0 

A b kinematics 

0.5 

q 2 migration 

0.4 

PID 

0.2 

Total 

+ 7.8 

- 8.2 


branching fractions can safely be considered in¬ 
dependent of the theoretical input at the current 
level of precision. 

From the ratio of yields and their determined 
efficiencies, the ratio of branching fractions of 
A° b -► pp~V^ to A° b -)• A+p~i 7 m in the selected q 2 
regions is 

g (AP+pp U M )q2 >15Ge vy c 2 

A+p-V fl ) q 2> 7Ge v/ c 2 

(1.00 ±0.04 ±0.08) x KT 2 , 


where the first uncertainty is statistical and 
the second is systematic. Using Eq. [I] with 
Rff = 0.68 ± 0.07, computed in Ref. |20| for 
the restricted q 2 regions, the measurement 



0.083 ±0.004 ±0.004, 


is obtained. The first uncertainty arises from 
the experimental measurement and the second is 
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due to the uncertainty in the LQCD prediction. 
Finally, using the world average \V c b\ = (39.5 ± 
0.8) x 1CU 3 measured using exclusive decays [id], 

| V u b | is measured as 

\V ub \ = (3.27 ± 0.15 ± 0.16 ± 0.06) x 10“ 3 , 

where the first uncertainty is due to the exper¬ 
imental measurement, the second arises from 
the uncertainty in the LQCD prediction and 
the third from the normalisation to \V c b\. As 
the measurement of V u b /1 already depends 
on LQCD calculations of the form factors it 
makes sense to normalise to the \V c b\ exclusive 
world average and not include the inclusive IQ, 
measurements. The experimental uncertainty is 
dominated by systematic effects, most of which 
will be improved with additional data by a reduc¬ 
tion of the statistical uncertainty of the control 
samples. 

The measured ratio of branching frac¬ 
tions can be extrapolated to the full q 2 re¬ 
gion using \Vcb\ and the form factor pre¬ 
dictions |20|, resulting in a measurement of 
B(A b —)■ pn~V^) = (4.1 ± 1.0) x 1CT 4 , where the 
uncertainty is dominated by knowledge of the 
form factors at low q 2 . 

The determination of \V u b\ from the mea¬ 
sured ratio of branching fractions depends on 
the size of a possible right-handed coupling |36|. 
This can clearly be seen in Fig. [4j which shows 
the experimental constraints on the left-handed 
coupling, \V^ b \ and the fractional right-handed 
coupling added to the SM, e#, for different mea¬ 
surements. The LHCb result presented here is 
compared to the world averages of the inclusive 
and exclusive measurements. Unlike the case for 
the pion in B° —> i t + £~V and B~ —> de¬ 

cays, the spin of the proton is non-zero, allowing 
an axial-vector current, which gives a different 
sensitivity to €r. The overlap of the bands from 
the previous measurements suggested a signifi¬ 
cant right-handed coupling but the inclusion of 
the LHCb \V u b\ measurement does not support 



Figure 4: Experimental constraints on the 
left-handed coupling, \V^ b \ and the fractional 
right-handed coupling, er. While the overlap 
of the 68% confidence level bands for the inclu- 
14 and exclusive [7j world averages of past 


sive 


measurements suggested a right handed coupling 
of significant magnitude, the inclusion of the LHCb 
\V ub \ measurement does not support this. 


that. 

In summary, a measurement of the ratio of 
V u b to \V c b\ is performed using the exclusive 
decay modes A b —> pp~v^ and A b —> A^^u^. 
Using a previously measured value of |, \V u b\ 
is determined precisely. The \V u b\ measurement 
is in agreement with the exclusively measured 
world average from Ref. |7|, but disagrees with 
the inclusive measurement 


14 


at a significance 
level of 3.5 standard deviations. The measure¬ 
ment will have a significant impact on the global 
fits to the parameters of the CKM matrix. 
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